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ABSTRACT 

Using A/'-body simulations we study the tidal evolution of initially disky dwarf galaxies orbiting a 
Milky Way-like host, a process known to lead to the formation of dwarf spheroidal galaxies. We focus 
on the effect of the orientation of the dwarf galaxy disk’s angular momentum with respect to the orbital 
one and find very strong dependence of the evolution on this parameter. We consider four different 
orientations: the exactly prograde, the exactly retrograde and two intermediate ones. Tidal evolution 
is strongest for the exactly prograde and weakest for the exactly retrograde orbit. In the prograde case 
the stellar component forms a strong bar and remains prolate until the end of the simulation, while 
its rotation is very quickly replaced by random motions of the stars. In the retrograde case the dwarf 
remains oblate, does not form a bar and loses rotation very slowly. In the two cases of intermediate 
orientation of the disk, the evolution is between the two extremes, suggesting a monotonic dependence 
on the inclination. We interpret the results in terms of the resonance between the angular velocity 
of the stars in the dwarf and its orbital motion by comparing the measurements from simulations 
to semi-analytic predictions. We conclude that resonant effects are the most important mechanism 
underlying the tidal evolution of disky dwarf galaxies. 

Subject headings: galaxies: dwarf — galaxies: fundamental parameters — galaxies: kinematics and 
dynamics — galaxies: structure — Local Group 


1. INTRODUCTION 

The formation of dwarf spheroidal (dSph) galaxies in 
the Local Group remains an open question but one of the 
most promising scenarios for their origin is via the tidal 
interaction of their disky progenitors with more massive 
hosts like the Milky Way. The scenario, proposed by 
Mayer et al. (2001) explains the morphology-density re¬ 
lation observed among the dwarfs of the Local Group 
and accounts for the non-sphericity of the dSph objects. 

The efficiency of the mechanism and its observational 
predictions have been investigated in detail by Klimen- 
towski et al. (2009), Kazantzidis et al. (2011) and Lokas 
et al. (2011, 2012). These studies explored the depen¬ 
dence of the process on a large number of orbital and 
structural parameters of the dwarf. The general picture 
that emerged from these studies is that a disky dwarf 
progenitor, once accreted by a massive host, undergoes 
strong tidal stirring and mass stripping if the orbit is 
tight enough. 

Typically, at the first pericenter passage, the disk 
transforms into a tidally induced bar (for a detailed de¬ 
scription of the properties of such a bar see Lokas et 
al. 2014a). The bar becomes thicker and shorter in time 
leading in the end to the formation of a spheroidal stellar 
component. The morphological transformation is accom¬ 
panied by strong changes in the kinematics as quantified 
by the amount of ordered to random motion. The lat¬ 
ter starts to dominate at some point and at the end the 
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galaxy is pressure supported. 

Among the parameters expected to have a strong im¬ 
pact on the evolution is the inclination between the an¬ 
gular momentum of the dwarf galaxy disk and its orbital 
angular momentum. However, in the studies mentioned 
above only a narrow range of inclinations was studied 
in detail, namely those with values i = 0°, 45° and 
90°. Because of this range, and the way the proper¬ 
ties of the dwarf were measured, no clear evidence for 
the dependence on this parameter was found. On the 
other hand, the difference between the prograde and ret¬ 
rograde galaxy encounters has been recognized for a long 
time and known to lead to very different outcomes (e.g. 
Holmberg 1941; Henon 1970; Kozlov et al. 1972; Toomre 
& Toomre 1972; Keenan & Innanen 1975). The issue has 
been recently addressed again by D’Onghia et al. (2009, 
2010) using an improved version of the impulse approx¬ 
imation applied to rotating systems. Although this ap¬ 
proximation is not directly applicable to our simulations, 
we attempt a comparison between these results and the 
numerical ones. 

In this paper we aim at clarifying the issue of the 
dependence of the results of tidal encounters between 
dwarfs and their hosts on the inclination of the dwarf’s 
disk. For this purpose we performed four simulations 
of tidal evolution of a dwarf galaxy orbiting a Milky 
Way-like host with disk inclinations i = 0°, 90°, 180° 
and 270°. The angles of 0° and 180° correspond to ex¬ 
actly prograde and exactly retrograde orientations of the 
dwarf’s disk. We also measured the properties of the 
dwarf galaxy in a different way that enables clear com¬ 
parisons between different runs. Preliminary results of 
this study, using lower resolution simulations, were dis¬ 
cussed in Lokas & Semczuk (2014). 

The paper is organized as follows. In section 2 we 
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TABLE 1 

Initial conditions for the simulations 


Simulation 

Ax 

Ay 

Az 

Inclination 

(deg) 

Line color/type 

10 

0.0 

0.0 

1.0 

0 

red/solid 

190 

0.0 

-1.0 

0.0 

90 

blue / short-dashed 

1180 

0.0 

0.0 

-1.0 

180 

green/dotted 

1270 

0.0 

1.0 

0.0 

270 

cyan/long-dashed 


present the simulations used in this study. In section 3 
we discuss the properties of the dwarf galaxies as they 
evolve in time, focusing on their kinematics, morphology 
and density profiles. Section 4 compares the results of the 
simulations to the predictions of semi-analytic models. 
The conclusions follow in section 5. 

2. THE SIMULATIONS 

The simulations used in this study were similar to those 
described in detail in Lokas et al. (2014a). Here we there¬ 
fore provide only a short summary. The initial condi¬ 
tions for the simulations consisted of A/'-body realizations 
of two galaxies: the Milky Way-like host and the dwarf 
galaxy, generated via procedures described in Widrow & 
Dubinski (2005) and Widrow et al. (2008). Both galaxies 
contained exponential disks embedded in NFW (Navarro 
et al. 1997) dark matter haloes, each made of 10 6 parti¬ 
cles (4 x 10 6 total). We note that the results presented 
here differ only slightly (are less noisy) from those of 
lower resolution simulations in Lokas & Semczuk (2014) 
where a smaller number of particles was used (by a fac¬ 
tor of five). We are therefore confident that our present 
resolution is sufficient to grasp all the essential features 
of the evolution. 

The dwarf galaxy model had a dark halo of mass Mh = 
10 9 M 0 and concentration c = 20. Its disk had a mass 
Md = 2 x 10 7 M 0 , an exponential scale-length Rd = 0.41 
kpc and thickness Zd/Rd = 0.2. The model is stable 
against formation of the bar in isolation for the time 
scales of interest here. The host galaxy was similar to 
the model MWb of Widrow & Dubinski (2005). It had 
a dark matter halo of mass Mh = 7.7 x 10 11 M 0 and 
concentration c = 27. The disk of the host had a mass 
Md = 3.4 x 10 10 M 0 , the scale-length R& = 2.82 kpc 
and thickness z d = 0.44 kpc. The disk was also stable 
against bar formation to avoid strong variations of the 
host potential in time. The disk of the Milky Way was 
coplanar with the orbit of the dwarf. Although this may 
seem contrary to observational constraints where most of 
satellite orbits are found to be polar (e.g. Pawlowski & 
Kroupa 2014), the choice was motivated by the necessity 
to avoid any additional variability which may be due to 
the passages through the plane of the Milky Way disk. 
However, we have performed an additional simulation to 
verify that for the orbits used here the evolution of the 
dwarf depends very weakly on the orientation of the orbit 
with respect to the Milky Way disk. 

The dwarf galaxy was initially placed at an apocen- 
ter of a typical, eccentric orbit around the Milky Way 
with apo- to pericenter distance ratio of r a po/r per i = 
120/25 kpc. The initial position was at the coordinates 
(X,Y,Z) = (—120,0,0) kpc of the simulation box and 


the velocity vector of the dwarf was toward the negative 
Y direction. We performed four simulations with differ¬ 
ent dwarf disk orientations with respect to the orbit: two 
coplanar with the orbit (prograde and retrograde) and 
two perpendicular to the orbit with angular momenta in 
the same and opposite direction to the dwarf’s orbital 
velocity. The different initial conditions, in particular 
the components of the unit angular momentum vector, 
are listed in Table [lj We will refer to the simulations by 
names indicating the initial inclination of the disk, 10, 
190, 1180 and 1270, where the inclination is measured as 
the rotation angle around the X axis of the simulation 
box. 

The evolution of the system in each simulation was 
followed for 10 Gyr using the GADGET-2 7V-body code 
(Springel et al. 2001; Springel 2005) with outputs saved 
every 0.05 Gyr. The adopted softening scales were = 
0.02 kpc and eh = 0.06 kpc for the disk and halo of the 
dwarf while cd = 0.05 kpc and en = 2 kpc for the disk 
and halo of the host, respectively. 


3. EVOLUTION OF THE DWARFS 

In this section we look at the inner properties of the 
dwarf galaxies as they are transformed by the tidal forces 
from the Milky Way. All measurements discussed below 
were made for stars (and dark matter particles) within 
the radius of 0.5 kpc from the center of the dwarf. 

3.1. Mass content 

We begin the analysis by measuring the mass inside 
this radius. Figure [T| compares the mass of stars (upper 
panel) and dark matter (lower panel) in the four simu¬ 
lations. As expected, the dark mass content decreases 
systematically, most significantly at pericenters that oc¬ 
cur at t = 1.2, 3.3, 5.5, 7.6, 9.7 Gyr from the start of 
the simulation and there is no dependence on the incli¬ 
nation because the halos were spherical and isotropic in 
all cases. 

On the other hand, for the stellar component (upper 
panel) there is a significant difference between dwarfs 
with varying initial inclination. However, the depen¬ 
dence may seem surprising because we expect the pro¬ 
grade disk to be more stripped, while the opposite is 
seen in the Figure: the stellar mass for the 10 simulation 
is even increased after the first pericenter passage, while 
it is decreased in the remaining cases. While the tidal 
stripping is indeed stronger in the outer parts, in the 
inner region we probe by this measurement, the stellar 
content increases as a result of the significant change in 
the structure of the stellar component due to the forma¬ 
tion of a tidally induced bar. 
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Fig. 1.— The evolution of the stellar (upper panel) and dark 
(lower panel) mass of the dwarf galaxy enclosed within the radius 
of 0.5 kpc. Different lines correspond to different initial inclination 
of the dwarf galaxy disk. 

3.2. Kinematics 

The main criteria usually applied in order to verify if a 
dwarf galaxy transformed from a disky object to a dSph 
are based on kinematics and shape of the stellar compo¬ 
nent. A dSph galaxy is supposed to be characterized by 
the dominance of random motions of the stars over the 
amount of rotation and its shape should be sufficiently 
close to spherical. In this and the following subsection 
we look in detail at the evolution of these quantities. 

In order to measure these properties, for each simula¬ 
tion output we determine the directions of the principal 
axes of the stellar component from stars within the ra¬ 
dius of 0.5 kpc using the inertia tensor and rotate these 
stars so that the new coordinate system is aligned with 
the principal axes (the x axis is along the longest, the 
y axis along the intermediate and the z axis along the 
shortest axis of the stellar component). We then intro¬ 
duce a standard spherical coordinate system such that 
< j) measures the angle around the z axis and 6 the angle 
from the z axis towards the xy plane. 

The kinematic properties were estimated using these 
coordinates. In the top panel of Figure [ 2 ] we plot the 
rotation around the shortest axis V = as a func¬ 
tion of time. We note that there is no other significant 
streaming motion along the other spherical coordinates 
or around the two other principal axes (see Lokas et al. 
2014b for a brief discussion of this issue). Clearly, the 
rotation is decreasing most strongly for the exactly pro¬ 
grade inclination of the disk ( 10 ), a significantly smaller 
decrease is seen for the perpendicular orientations (190 
and 1270) and the effect is the weakest for the retro¬ 
grade case (1180). Note that the decrease of rotation is 
not steadily monotonic even in the exactly prograde case 
(10). This is due to the tidal torques acting on the bar 
at pericenter passages that can speed up or slow down 
the bar depending on its particular orientation at this 



time [Gyr] 


Fig. 2.— The evolution of the mean rotation velocity (upper 
panel), the velocity dispersion (middle panel) and the ratio of the 
two (lower panel) as a function of time. 

moment (see Lokas et al. 2014a for details). 

The behavior of the velocity dispersion is exactly the 
opposite. In the middle panel of Figure[2]we show the ID 
velocity dispersion calculated as a = [(cr^H-cr^H-o - ^)/3] 1//2 . 
The increase of a at the first pericenter passage is 
strongest for the 10 case, intermediate for 190 and 1270 
and for 1180 a remains constant in time or even slightly 
decreases due to mass loss. The ratio V/a shown in the 
lower panel of Figure [2] decreases for all simulations, but 
reaches a value significantly below unity only for the pro¬ 
grade case. For the three remaining cases a substantial 
amount of rotation is retained, although the hierarchy of 
lower V/a for more prograde cases is preserved. 

In Figure [3] we show the evolution of the different veloc¬ 
ity dispersions a r (upper panel), oq (second panel) and 
CT 0 (third panel) as a function of time. In all cases the 
increase of a given dispersion at the first pericenter pas¬ 
sage and its level at later times is highest for the 10 case, 
intermediate for 190 and 1270 and non-existent for 1180. 
In addition, this increase is the most abrupt for the pro¬ 
grade dwarf, while for the intermediate inclinations 190 
and 1270 the increase occurs much more slowly in time 
and takes about half the orbital period between the first 
and second pericenter passage. Interestingly, significant 
increase is seen in all dispersions in spite of the fact that, 
due to the formation of the bar, one could expect the 
radial a r to increase much more significantly. 
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time [Gyr] 


Fig. 3.— The evolution of the velocity dispersions of the stars 
in spherical coordinates cr r (upper panel), ctq (second panel) and 
(7^ (third panel). The lower panel shows the evolution of the 
anisotropy parameter j3. 

We end the analysis of the kinematics by plotting in 
the lower panel of Figure [3] the value of the anisotropy 
parameter (3 = 1 — erf/(2of) where erf = erf + cr 2 + V£ is 
the tangential second moment including rotation. A sys¬ 
tematic dependence of this parameter on the inclination 
of the disk is also present: the stellar orbits are most ra¬ 
dial for the 10 case with (3 almost constantly at the level 
of 0.5, corresponding to mildly radial orbits of the stars, 
characteristic of the bar. For the intermediate cases 190 
and 1270 the (3 values stay between 0 and 0.4, while for 
the retrograde case 1180 f) remains negative due to the 
dominant presence of rotation. Interestingly, at the end 
of the evolution all three non-prograde cases 190, 1180 
and 1270 have almost isotropic orbits (/? = 0) although 
they reach this special value via different evolutionary 
paths. 



time [Gyr] 


Fig. 4.— The evolution of the shape of the stellar component 
in time. The panels from top show the axis ratios 6/a, c/a, the 
triaxiality parameter T and the bar mode A 2 . 

3.3. Shapes 

Figure [4| illustrates the evolution of the shape of the 
stellar component of the dwarfs in time. In the first and 
second panels from top we plot the axis ratios b/a (inter¬ 
mediate to longest) and c/a (shortest to longest). The 
thickening of the dwarf, as quantified by the increasing 
value of c/a is similar in all cases, although for the pro¬ 
grade case c/a decreases after the first pericenter due to 
the formation of the bar. Much more significant differ¬ 
ences are seen in the evolution of b/a. This value de¬ 
creases most strongly for the prograde 10 case signifying 
the prolate shape characteristic of the bar. The interme¬ 
diate inclinations 190 and 1270 lead to less prolate shapes, 
while 1180 remains disky for the whole evolution, as in¬ 
dicated by b/a remaining constantly close to unity. 

The shape can be also quantified in terms of the tri- 
axiality parameter T = [1 — (b/a) 2 ]/[l — (c/a) 2 ] which 
is shown in the third panel of Figure 01 The values of 
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Fig. 5.— Surface density distributions of the stars in the dwarfs 
at the last apocenter (t = 8.65 Gyr) seen along the intermediate (y) 
axis of the stellar component. The surface density measurements 
were normalized to the maximum value X max = 4.8 X 10 5 stars 
kpc -2 occurring for 1180. Contours are equally spaced in logE 
with AlogX = 0.05. 


the parameter T > 2/3 at all times confirm the prolate 
shape due to the bar in the case of 10. For the intermedi¬ 
ate cases 190 and 1270 we get 1/3 < T < 2/3 indicating a 
triaxial shape. For the retrograde case 1180 the parame¬ 
ter remains low, T < 0.2, at all times as is characteristic 
of a disk. 

The presence of a bar is usually detected by measur¬ 
ing the bar mode A 2 of the Fourier decomposition of the 
stars projected along the shortest axis of the stellar dis¬ 
tribution (see a more detailed discussion in Lokas et al. 
2014a). Usually, A 2 > 0.3 is considered as high enough 
to be interpreted as a bar. As we can see in the lower 
panel of Figure [4| this is always the case for simulation 10 
after the first pericenter passage and also for some signif¬ 
icant periods of time for the intermediate cases 190 and 
1270, which means that the bar also forms there, but it is 
much weaker. The exactly retrograde disk does not form 
a bar as its A 2 < 0.06 at all times. Slight temporary 
increases of this value are due to stretching of the dwarf 
at the pericenters. 

These measurements are confirmed by the maps of the 
surface density distribution of the stars in the four dwarfs 
plotted in Figured The distributions are shown in pro¬ 
jection along the intermediate axis of the stellar com¬ 
ponent so in their most non-spherical appearance. The 
snapshots were selected for the time t = 8.65 Gyr af¬ 
ter the start of the simulations, corresponding to the last 
apocenter passage. In all cases, except for the retrograde 
one 1180, the remnant of the bar formed after the first 
pericenter passage is still visible in the inner parts. In ad¬ 
dition, the distribution of the stars in the prograde case 
10 is much less diffuse in the outer parts of the maps. 
In spite of the fact that this dwarf evolves most strongly, 



. . . . 

0.1 1 10 

r [kpc] 


Fig. 6.— Comparison of stellar density profiles of the dwarfs at 
the fourth apocenter (6.6 Gyr after the start of the simulations, 
colored lines). For completeness we also show the dark matter 
profile (black line) which is similar for all simulations. 

its neighborhood is not uniformly filled with stripped de¬ 
bris because the lost stars form well-defined, narrow tidal 
streams. 


3.4. Density profiles 

The dependence of the properties of the dwarfs on the 
initial inclination of the disk also manifests itself in the 
evolved density profiles. In Figure [6] we show examples of 
the stellar density profiles (different colors) and the dark 
matter density profile (black line) for different simula¬ 
tions considered here, measured at the fourth apocenter. 
The transition from the bound component to the tidal 
tails is visible as the break in the slope of the density 
profiles. This transition is however only well defined for 
the exactly prograde case 10 (and the dark matter pro¬ 
file). In this case the transition from the steeper to the 
shallower profile (where the slope is around r -4 ) occurs 
at around 3 kpc. For other simulations no such clear 
break radius is seen. 

As discussed in Lokas et al. (2013) using similar (but 
only mildly prograde) simulation setups, the break radii 
can be interpreted as the tidal radii. In this case, the 
dependence on the orbit of the star within the satellite 
is expected (Keenan & Innanen 1975; Read et al. 2006) 
and we will attempt a detailed comparison in a follow-up 
paper. Since the stellar density profiles of our simulated 
dwarfs do not show clear signatures of the break radius 
(except for the exactly prograde case) here we propose 
a comparison in terms of density. The stellar profiles 
shown in Figure [6] demonstrate clear hierarchy: at the 
outer radii (larger than 1 kpc) the stellar density profile 
of 1180 (green line) is above all the other profiles, the one 
of 10 (red lines) is the lowest, and the ones of 190 and 1270 
fall exactly on top of each other and between the other 
two. This means that we can quantify the amount of 
tidal stripping in these different cases by measuring the 
density of the stars, rather than the break radius. 

To do so, we calculated the mean density of stars in the 
shells of radii 2 kpc < r < 4 kpc as a function of time for 
different simulations. The results are shown in the upper 
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Fig. 7.— Upper panel: the density of stars in the shell of radii 2 
kpc < r < 4 kpc as a function of time. Lower panel: the radius at 
which the density of stars drops below 10 4 MQ/kpc 3 as a function 
of time. 

panel of Figure [71 In the lower panel we plot the radius 
rd at which the density of stars drops below lO 4 M 0 /kpc 3 
as a function of time. In both plots there is a clear sys¬ 
tematic difference between the measurements for differ¬ 
ent initial disk orientations: the stars are stripped more 
effectively on prograde orbits as demonstrated by the 
lower outer densities and smaller radii where the density 
drops to a fixed value. Although the values for the mea¬ 
surements were chosen in an arbitrary way, we expect 
the results to be similar if these parameters are slightly 
modified. They clearly confirm that the amount of tidal 
stripping depends very strongly on the initial inclination 
of the dwarf’s disk. 


4. COMPARISON WITH SEMI-ANALYTIC PREDICTIONS 

To describe the encounter between a disky dwarf 
galaxy and the Milky Way we use the impulse approx¬ 
imation as discussed in Binney & Tremaine (1987) and 
D’Onghia et al. (2010). According to this approximation 
the kth component of the acceleration of each dwarf’s 
star with respect to the center of mass is given by 


Vk = ~ 


E 


3 


dxkdxj 


X j , 

£C —0 


( 1 ) 


where is the Milky Way’s (i.e. perturber’s) potential 
and Xi are Cartesian coordinates. We integrate equation 
m over a finite time period to obtain velocity increments 
that can be compared with increments measured from 
simulations 

At 

Avfc = J i>k d t. (2) 

o 

We work in a frame centered on the center of the 
dwarf’s mass, the dwarf’s disk lies in xy plane and po¬ 
sition of the Milky Way at the pericenter is (b x ,0,b z ). 
In this frame the trajectory of each star originating from 


the dwarf is 

x = [r cos (Clt + </>o), r sin (Qt + 0o), 0], (3) 

where r is the radius, Q is the angular velocity and (j> o is 
the initial azimuthal angle of the star. The trajectory of 
Milky Way is given by 

X = Vt + X o, (4) 

where V and Xq are constant vectors, fitted to mimic 
the perturber’s trajectory from simulations as a straight 
line during a given time period. Note that V and X 0 
are different for each of our simulations as they depend 
on the inclination i. This dependence can be found by 
rotating the trajectory for the prograde case with matrix 
A defined by Euler angles, in order to obtain trajectories 
in other cases. One of the Euler angles is the inclination 
i and other two depend on the perturber’s orbit. 

To approximate the gravitational potential of the 
Milky Way we sum the potential from its stars and the 
dark matter halo. The first part is represented as a point- 
mass potential, while the second is given by the NFW 
profile 

GM d GMn ln(l + c\x — X\/r v ) 

f = -W^x]- s -F^xj-• <5) 

where Md is the mass of the Milky Way disk, Mr is the 
virial mass of its halo, r v is the virial radius, c is the 
concentration parameter and g = [ln(l + c) — c/( 1 + c )] -1 
(see Lokas & Mamon 2001). Substituting equations (pQ), 
©, © and m into m we obtain formulae which can 
be numerically integrated to get velocity increments. 

In order to calculate the velocity increments described 
above we estimate the angular velocity in the dwarf as 
Q = \v\/r. We calculate integrals of equation (| 2 j) over a 
small period of time, so that the assumptions concerning 
the trajectories are valid. In our simulations the outputs 
were saved every 0.05 Gyr and we choose to integrate 
over this time to compare our predictions with velocity 
increments occurring in the simulations when the tidal 
force is the strongest, i.e. between the output preceding 
the first pericenter and the one as close as possible to 
this pericenter. We check how the distribution of incre¬ 
ments changes with the distance from the center of the 
dwarf. The results for all simulations are summarized in 
Figure [HI 

The red points in Figure [ 8 ] represent the values pre¬ 
dicted by the impulse approximation, while the blue 
points correspond to the values measured from the sim¬ 
ulations. At radii smaller than 2 kpc the semi-analytical 
predictions do not reproduce the distribution of velocity 
increments because at these radii the velocity changes 
are dominated by the gravitational potential of the dwarf 
which was not included in the predictions. However, for 
radii greater than 2 kpc the agreement between the full 
simulations and semi-analytic predictions for this short 
time period is very good. In particular, in the upper left 
panel of Figure [ 8 ] for simulation 10 we find two very well- 
defined branches corresponding to stars on different sides 
of the dwarf. Branches from simulations are not exactly 
symmetric with respect to zero while the branches from 
theoretical predictions are. The difference is due to the 
fact that the analytic predictions only take into account 
the lowest order terms. 
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Fig. 8.— The components of velocity increments along the x, y and 2 : axis for simulations with different initial inclination of the dwarf’s 
disk. The blue dots indicate the values measured for a subsample of stars in the simulation. Red dots show the corresponding semi-analytic 
predictions. 


For simulations 1270 the increments are almost identi¬ 
cal as for 190, as expected due to symmetry of the two 
configurations with respect to the orbital plane. In some 
of the panels in Figure [8] the velocity increments for radii 
larger than 2 kpc are approximately zero. However, the 
negligible values are consistently obtained both from the 
simulations and the semi-analytic calculations. 

We further illustrate these results in Figure [9] where we 
plot the surface density maps of the stellar component 
and the velocity vectors for a random sample of a hun¬ 
dred stars at radii larger than 2 kpc. The plots show the 


dwarfs at the first pericenter passage, i.e. after they have 
been affected by a tidal impulse from the Milky Way for 
the first time. The coordinate system is as defined above, 
with xy coordinates in the plane of the dwarf’s disk. 

The comparison of the upper left panel of Figure [9j 
corresponding to simulation 10, to the other three con¬ 
firms that for the prograde encounter the effect of the 
tidal force in the strongest: the dwarf galaxy disk is al¬ 
ready strongly distorted toward a bar-like shape and two 
tidal arms are formed. There are significant increments 
of velocity along the x axis. Small increments of veloci- 
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Fig. 9.— Surface density distributions of the stars in the dwarfs at the first pericenter (t = 1.15 Gyr) in projection onto the initial disk 
plane. Arrows indicate the velocities of the randomly selected subsample of 100 stars outside the radius of 2 kpc. 


ties are also visible along the y axis in the right panels of 
Figure [9] corresponding to simulations 190 and 1270. The 
effect of the tidal force is weakest for the retrograde case 
(lower left panel of Figure [9]) where the dwarf’s initial 
disk is affected very little. 

As discussed by D’Onghia et al. (2010) strongest tidal 
interactions occur when the intrinsic angular velocities of 
stars in the dwarf’s disk are comparable to the angular 
velocity of the satellite on its orbit. This condition can 
be written as 

^disk — ^orb (0) 

and we can define the resonance parameter 


that should be of the order of unity for the strongest, 
resonant response. 

To demonstrate the resonant nature of the tidal effects 
in our simulations we measured a(r) for each simulation 
output and found radii from the center of the dwarf at 


which a = 1. The time dependence of this radius for 
simulation 10 is shown in Figure [lOl We can see that 
the variability of this radius reflects the varying orbital 
velocity fl 0 rb of the dwarf resulting in rather large val¬ 
ues at apocenters and much smaller ones at pericenters. 
The slow variation over time scales much larger than the 
orbital period is caused by the mass loss and decreas¬ 
ing ft. Whenever the dwarf reaches the pericenter of its 
orbit around the Milky Way, this characteristic radius 
drops below 2 kpc, and even down to 1 kpc at the later 
pericenters. Note that these radii are of the order of 2- 
3 half-light radii of the stellar component which means 
that a significant fraction of stars is affected. Comparing 
with Figure [8] showing the velocities at the first pericen¬ 
ter we confirm that this radius (equal to 1.7 kpc at this 
time) is exactly where the tidal effects start to prevail 
over the dwarf’s potential. 

5. CONCLUSIONS 

In this work we extended previous studies of the effi¬ 
ciency of the tidal stirring mechanism to include the de- 
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Fig. 10.— The evolution of the radius at which the resonance 
parameter a = 1 in time. Only results for simulation 10 are shown 
because the measurements yield similar values for the remaining 
simulations. 


pendence on the initial inclination of the dwarf galaxy 
disk with respect to its orbit around the Milky Way. 
Our simulation setups involved a dwarf galaxy placed 
on a typical, eccentric orbit around a Milky Way-like 
host and its evolution was followed for 10 Gyr. We con¬ 
sidered four configurations, an exactly prograde, an ex¬ 
actly retrograde and two intermediate orientations of the 
disk. We found the efficiency of tidal stirring to be very 
strongly dependent on this inclination. 

The effect of the tidal interaction turns out to be 
the strongest for the exactly prograde orientation of the 
dwarf’s disk (10). In this case the disk transforms into a 
strong bar (A 2 = 0.4) at the first pericenter passage and 
a remnant bar is retained until the end of the evolution. 
Although the bar becomes weaker with time, the shape of 
the stellar component is consistently prolate, but tending 
to a spherical toward the end of the evolution. This mor¬ 
phological transformation is mirrored in the kinematics 
by the gradual decrease of the dwarf’s rotation veloc¬ 
ity. Already at the first pericenter passage the rotation 
drops significantly and the velocity dispersion increases, 
mostly in the radial direction due to the formation of the 
bar. During subsequent pericenter passages the rotation 
is further decreased down to rather low values at the end 
of the simulation where V/a = 0.3. Thus the stream¬ 
ing motions are almost completely replaced by random 
motions of the stars. Interestingly, the radial velocity 
dispersion dominates the whole time, which manifests it¬ 
self in the anisotropy parameter close to f3 = 0.5 even at 
the end. 

In the two cases of perpendicular orientations of the 
dwarf’s disk with respect to the orbit (190 and 1270) the 
evolution is to some extent similar to the 10 case. In these 
configurations the bar also forms after the first pericen¬ 
ter passage, but more slowly (over a time scale of about 
1 Gyr corresponding to half the orbital period) and is 
significantly weaker (A 2 = 0.3), also in the subsequent 
evolution. The overall shape of the stellar component can 
be characterized more as triaxial than decidedly prolate. 
The transition from the streaming to random motions 
of the stars also happens less efficiently with V/a only 
slightly below unity at the final outputs. The anisotropy 
parameter is close to zero due to the contribution of the 
still significant rotation. 


In the exactly retrograde case (1180) no strong evolu¬ 
tion is present: the dwarf’s stellar component does not 
form a bar and remains disky. The only signatures of 
tidal evolution in this case are the mass loss (similar as 
in other cases, mostly in dark matter), small decrease of 
the rotation velocity, slight evolution of the anisotropy 
parameter from negative toward isotropic and a non- 
negligible thickening of the disk. The difference between 
this and the other cases is also visible in the stellar den¬ 
sity profiles which are less affected and do not show any 
clear transition from the bound component to the tidal 
tails. 

We have interpreted these changes in the context of 
the resonant stripping mechanism recently discussed by 
D’Onghia et al. (2009, 2010). In particular, we calcu¬ 
lated the velocity increments the dwarf’s stars should ex¬ 
perience at the first pericenter and compared them with 
the direct measurements from simulations. We find a 
very good agreement between the two, confirming the 
interpretation that the evolution we see in the full N- 
body treatment is indeed due to the orientation of the 
dwarf’s disk. The resonant nature of the phenomenon is 
further confirmed by the behavior of the ratio between 
the angular velocity of the stars in the dwarf and the an¬ 
gular velocity of its orbital motion. This ratio turns out 
to be of the order of unity only near the pericenters and 
this is indeed when the tidal effects are the strongest. 

The results presented here suggest that the most im¬ 
portant mechanism underlying the tidal evolution of 
disky dwarfs orbiting a bigger galaxy is indeed of res¬ 
onant nature. We propose to refer to the processes of 
morphological and dynamical evolution of the dwarfs we 
described as ‘resonant stirring’ in analogy to the ‘res¬ 
onant stripping’ mechanism found by D’Onghia et al. 
(2009, 2010) to increase the mass loss in similar configu¬ 
rations. As discussed by D’Onghia et al. (2010), the res¬ 
onance is broad, hence the name ‘quasi-resonant stirring’ 
would be more appropriate. In physical terms, this reso¬ 
nance can be traced to the fact that the stars with a ~ 1 
remain for an extended period of time on the line join¬ 
ing the dwarf galaxy and the perturber. For these stars 
the tidal force (which is strongest along this line) has the 
longest time to operate which results in the largest veloc¬ 
ity increments and the largest stirring. In the context of 
the tidal radius calculations, the difference between the 
prograde and retrograde cases comes from the change of 
sign of the Coriolis force. The relation between the two 
approaches remains to be investigated and we plan to 
address this issue in our future work. 
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